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Abstract 

The attempt to describe the bell-shape dependence of the critical temperature of 
high-T c superconductors on charge carriers density is made. Its linear increase in the 
region of small densities (underdoped regime) is proposed to explain by the role of 
the order parameter phase 2D fluctuations which become less at this density growth. 
The critical temperature suppression in the region of large carrier densities (overdoped 
regime) is connected with the appearance (because of doping) of the essential damping 
of long-wave bosons which in the frame of the model proposed define the mechanism 
of indirect inter-fermion attraction. 

PACS: 67.20. +k, 74.20.-z, 74.25.-q, 74.72.-h 
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1 Introduction 



In spite of non-falling down the intensity of study of high-T c superconductors (HTSCs) the 
comprehension of many regularities which define their physical behaviour is not still achieved. 
First of all, it is put down to the HTSC normal state properties, the description of which in 
the framework of standard Fermi-liquid theory proves to be impossible if the carrier concen- 
tration rif in corresponding samples is like that one that the resulting critical temperature 
T c in them is less than some characteristic for every compound optimal value T c mai '. The 
latter in all (practically, without exclusion) HTSC copper oxides appears because, as it is 
well-known (see, for example, the rewiev [[I]]), the conductivity (consequently - superconduc- 
tivity) these, initially quasi-2D antiferromagnetic insulators, results from their doping by 
donor (Nd, Pr) or acceptor (Sr, Ba, O) ions. 

Then the temperature T c being formed at increasing of carrier density (or, what is the 
same, Fermi-energy) in a system also grows up. This growth, however, is rather quickly 
stopped and after some (not large upon nf values) part of "saturation" the drop to zero of 
the function T c (rif) follows, which thus acquires the bell-shape form. As a result, the HTSC 
compounds with relatively large concentration of itinerant carriers (so called, overdoped 
regime) become non-superconducting metals the behaviour of which to a certain degree can 
be carried out on the base of conventional theory of Fermi-liquid ||. The corresponding 
carrier concentration values for underdoped regime are such ones that the T c (rif) increase is 
observed for rif from ~ 0.04 — 0.08 carrier per Cu02 layer cell (i.e. from the threshold of the 
insulator-metal transition) to T™ ax ~ T c (0.15 — 0.18); the temperature T c (rif) becomes zero 
when rif ~ 0.25 — 0.28 || |], If the initial part of the function T c (rif), where dT c (n /) / dn f ~ 
const > H 0], can be (at least qualitatively) interpreted (see Refs.|8|, |9], [L(| and the review 
|TT| ) proceeding from the theory of the crossover from the Bose-Einstein condensation of 
separate composite bosons (local pairs) and their superfluidity to the superconductivity of 
the BCS type, then the reason of T c (rif) suppression in the BCS carrier concentration region 
remains, in fact, unknown yet. 

From the very begining the above-mentioned behaviour of T c (rif) was ascribed to the 
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strong electron-electron correlations (see, for example, the book ||T2|) and to the filling of 
usually narrow Hubbard subbands by itinerant carriers. This question was intensively inves- 



tigated by Prof. I.V. Stasyuk and his collaborators (13], [L4|, [15], |16 |. 

However a little bit later it becomes evident that carrier concentrations corresponding 
to a disappearance of superconductivity are so small that such a filling of the whole valent 
(in fact, conduction) band in HTSCs is impossible. The alternative version of the HTSC 
effect disappearance most consistently considered in Refs. fl7], |i~8H (see also fl], U) consists 



in assumption that bare Fermi-level in HTSCs proves to be in the vicinity of the extended 
saddle point of the electronic spectrum. According with this poin of view at rif growth 
the Fermi-energy €f — k^/2m (kp - Fermi momentum, m - effective mass) passes the van 
Hove singularity in the conduction band density of states and, hence, T c (rif) reveals its 
maximum. If it would be so the function T c (n/), being similar to BCS one, would have the 
exponential dependence on nf far from T c max what, as it is seen from experiments and was 
already said, is not confirmed by observations. On the other hand, the very appearance 
of the van Hove spectrum singularity requires fine-tuning (for example, because of rather 
definite ratio between the hole next-near and nearest hopping constants) what seems can be 
hardly achieved even in one compound not speaking of many ones. 

At the same time less attention was paid to the idea that the weakness of supercon- 
ductivity can be conditioned by "feedback" effect, or the changes in the spectrum u;(k) of 
intermediate bosons the exchange of which results in the fermion pairing. These bosons in 
HTSCs can be attributed to phonons (similarly to BCS-Eliashberg model) the role of which 



is actively advocated by Ginzburg (see, for example, his rewiev [0), magnons, or spin fluc- 



tuations, 0, ^TJ, quadrupole cM-excitons (Gaididei-Loktev- Weber mechanism (see Ref. 



plasmons etc. It must be, however, noted that if phonon spectrum does not (or very 
weakly) depend on doping then long-wave excitations of magnetic subsystem (cM-excitons, as 
excitations over a magnetic background, can be also considered here) are strongly suppressed 
because of the long-range magnetic order destruction in metallic phase of HTSCs [[I], p^ . 

Namely spin waves are subjected the most appreciable doping effect; their low-frequency 
region (for the wave vectors less than some characteristic value k m i n ) acquires diffusion 



form, or becomes overdamped. In the insulating phase of HTSCs as it is shown in Ref. |[23|| . 
kmin ~ 7i/; in their mettalic phase the magnon damping becomes even more strong and 
~ 71/ [||, U|, what is completely in line with degradation of the correlation 
magnetic length £ TOa „ ~ nT 1 ^ 2 measured in HTSCs p5|. In that way the long-wave damping 
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7(k) increases up to so high values (7(k) >> oo(k)) that corresponding intermediate bosons 
(most probably, spin fluctuations) can not participate in interaction transfer, "being out". 

The solution of the self-consistent magneto-electronic problem as a whole (i.e. the explicit 
solution of the equation of superconductivity similar to the Eliashberg equation) is scarcely 
possible now. Therefore below the attempt is made to consider the simplest model with 
indirect iner-fermion attraction what is provided by intermediate massive bosons with long- 
wave cut which is proportional to kp. The carriers to be supposed appear in a system 
due to doping what corresponds to generally accepted scenario of metallization of copper 
oxides. For the sake of simplicity we shall omit any other dampings (in particular, the carrier 
damping because of disorder the effect of which was analysed, for example, in Ref. |[26| ) . 

Following to Refs.f^, |Tl| we shall also suppose that in a model 2D system (in fact, all 



HTSCs can be with a good accuracy put down to this kind) the superconducting condensate 
is formed in a way principally different from the ordinary one. In such a case one has to differ 
the order parameter formation temperature from the real (observable) critical temperature 
in a sense that the absolute value of the former does not become zero at and above T c . 

2 Model and main equations 

The model Hamiltonian density of the electron-phonon system can be written in the well- 
known form: 

H(x) = -iPl(x) + f^j Vv(x) + H ph (<p(x)) 

+g ph ^t(x)^ a (x)ip(x), x = r,t, (1) 

where H p h is the Hamiltonian of free phonons which will be described more precisely below; 
ip a {x) and ip(x) are the fermion and boson field operators, respectively; m is the effective 



mass of the fermi-particles, a ="[*, J, - their spin variable, and g p h is the electron-phonon 



density in a system; we put % = k B = 1. 

As it was mentioned in Introduction, we shall model boson-exchange interaction taking 
into account the dependent on doping saturation of long-wave bosons. Under above men- 
tioned assumptions about gradual non-participation of the part of the bosons in an attracting 
interaction formation the simplest way to describe that can be achieved by making use of 
the free-phonon propagator in the form: 



where u;(k) is the boson dispersion law, and 6{k) is the step function. As it was pointed out 
above k mi „ in (0) is some characteristic wave vector which devides the region of overdamped 
(k < k min ) and long-lived (k > k min ) intermediate bosons (here: phonons by definition). 
Although magnetic correlation length measurements result in k m i n = 2kp we shall adopt 
a more general (or soft) relation supposing that k min = akp = ay/2mep where a is some 
free-parameter. 

It is very important that the Hamiltonian (|I]) is invariant with respect to the global 
symmetry transformations 



which in 2D case (unlike to 3D one) remain unbroken and phase transition is here accompa- 
nied by change in the correlation function behaviour only. 

The T—rif phase diagram of a system can be calculated by using the Hubbard-Stratonovich 
method, generalized on the case of non-local (indirect) interaction (so called auxiliary bilocal 
field method). For the finding of the grand partition function Z it is useful to pass to Nambu 
spinors: \Ed — (V'jjV'l) an d its conjugated one. After that performing the integration over 
bosonic fields it is easy to obtain the Lagrangian 



coupling constant; fi in ([!]) is the chemical potential of fermions, which fixes their average 




(2) 




(3) 



L = tf(x)[-d T + ( — + fi)r z ]^(x) 




(4) 
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of a system, where the integration over repeated indices is supposed. The kernel K in (B) 
will be defined below. 

Let introduce the pairing order parameter 

<fi(xi, Vi) = K(x u yi ; x 2 , S/a)*^)*^)?"* 

= r + 0(xi,yi) +r_^*(zi,yi), (5) 

where r + = |(t x + iT y ), r_ = ^(t^ — ir y ) (and r 2 in are the Pauli matrices. 
Then adding to the Lagrangian L a zero term 

^[4>(xi, Vi) - K(xi,yi, x[, y l 1 )^f(x' 1 )^f i (y' 1 )T z ]K~ 1 (x 1 , y x \ x 2 , y 2 ) [<f>(x 2 , y 2 ) 

-K(x 2 , y 2 ; x' 2 , y' 2 )^{x' 2 )^{y' 2 )T z ] 

with the purpose to cancel the four-fermion interaction one comes to the expression: 

V 2 

L(x 1 ,y 1 ;x 2 ,y 2 ) = W{x{)[-d T + (- h /x)r 2 - T + (f>(x u yi) 

2m 

-r-(f)*(xi, 2/i)]*(yi) + -</>0i, yi) yi; x 2 , y 2 )<fi(x 2 , y 2 ) (6) 

for the Lagrangian needed. The Fourier transformation of K can be written as 

d 3 Pd 3 p 1 d 3 p 2 ,xi+yi x 2 + y 2 



K(x l ,y 1 ;x 2 ,y 2 ) = j —-- K P {p 1 ;p 2 )exp[-iP( 



(27T) 

-zpi(xi - yi) - ip 2 (x 2 - y 2 )}, 

(pi = (p i , cjj) where i — 1,2 and P = (P, u;) are the relative and the centre of mass momenta, 
respectively). Supposing now that Kp(pi;p 2 ) is P- independent we pass to the standard 
kernel form 

K(pi,p 2 ) = g 2 ph D{p l -p 2 ), (7) 

which corresponds to the indirect inter-fermion interaction. 
The partition function can be written then as 

Z = J V^V^V(f)V(f)* exp -f3 J L(^\ty,(f)*,(f))dxdy 

= J P0£>0*exp(-/ft)[0]), p = 1/T, 



In the last expression Q[Q] is the thermodynamic potential which in the "leading order" on 
g p h has the form 



pn[g] 



-Tr 



Ln^- 1 + i(^jr-V) 



in which 



»-i 



(8) 



(9) 



is the full fermion Green function. After the direct minimization of the potential (||) it is 
easy to obtain the equation for auxiliary 0-field (or order parameter): 

(Pkduj 



5n/5(P = 20 - tr J -F^-K(p; k, u)Q(k, u)r x = 0. 



(2k) 

Using ([T^) one can readily arrive to the well-known Cornwell-Jackiw-Tombou- 
lis formula for the effective action in the one- loop approximation P^ : 

pn(G) = -TrLn^ 1 + ^TrOKQ, 



(10) 



or, taking into account (|9]) (or (|T0"D), 

f3Q(G) = -Tr 



i^+^o 1 -!; 



To investigate the possibility of the condensate formation in a 2D system it is convenient 
in accordance with || to pass to a modulus-phase parametrization of the order parameter 
(cp. @, where 9 = const): 



<f>(x, y) = p( x , y)exp[-i(0(x) + 0(y))/2] 



(12) 



with simultaneous Nambu spinor transformation 



\f ''(x) = x (x)exp[i0(x)T z /2], 



(13) 



corersponding to "separation" of the bare fermions on their neutral x( x ) an d charge 8(x) 
parts (fermi- and bose-ones, respectively). 

In the approximation that p(x,y) — p — const (see Ref . [p8fl ) and spatial ^-fluctuations 
are small one can (using ( |T2| ) and (|T3|)) obtain the next expressions for Q and Q which are 
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defined in (ED and (HIT): 



g- 1 



-r z T 9 + 



w 2 



— I 



\ 2 e wv 1 



2m 



m 



2m , 

G-\ P )-v{de) (14) 



and f2 = flkin{p, V0) + ^pot(p) with fl po t(p) which is defined by ( |TTD at V6* = and 
{3n km (p, V0) = Tr[GS - G S + ^GSGS - ^G ZG Z 

+ T x ^ipG{GE + GEGE)) = ^ dr J d 2 vJ(p,T, p(p,T))(V9) 2 , (15) 
where the effective neutral fermion stiffness 



1 

87 



p 2 + p 2 + p + 2T\n 1 + exp( 



T 

An 







AT 2 d(p 2 /AT 2 ) 



dx- 



Vp 2 + p\ 

T } 
x + p/2T 



(16) 



-m/2T cosh 2 ^x 2 + p 2 /AT 2 
was introduced. 

The evident analogy with XY-model (two-component order parameter in 2D space) gives 
the equation for the temperature T BKT of the Berezinskii-Kosterlitz-Thouless phase transi- 
tion in a system, namely (see Chapter 15 in the book |p9|| ): 



7T 



— J(p,Tbkt, p(p,Tbkt)) — Tbkt- 



(17) 



(Recall that the temperature Tbkt plays the role of critical one in 2D metals.) 

The parameters p and p in ( |T7]) being dependent of T are still unknown; therefore it is 
necessary to obtain the equations which connect them with carriers density nj. The first 
one follows from (|T0D for p ^ 0: 

d 2 k K{u m ) 



1 = T £ 



(2n) 2 u 2 m + e(k)+p 2 
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where u n = (2n + l)irT is the Matsubara fermion frequencies, £(k) = k 2 /2m — p and it was 
used the Einstein model for phonon dispersion: ui(k) = uq. The dependence on parameter 
a in Eg. (pop is preserved through the kernel K (see (7) and (2)). 

The second one is defined by the condition V~ l dfl pot (p)/dp = —nf (V is the volume of 
a system) what results in the well-known number equation: 



p 2 + p 2 + p + 2T In 



1 + exp 



2e F , (19) 



T 

where the equality e F = 7mf/m was used; it is correct for free 2D fermions with a quadratic 
dispersion law. 

Thus, we have obtained a self-consistent set of the Eqs.(17)-(19) needed to investigate 
the phase diagram of a 2D metal with arbitrary carrier density. The latter is the parameter 
defined such a metal superconducting properties. 



3 Phase diagram of a system 

As it can be seen from the previous Section (see also Refs.||, |ll|]) there exist two characteristic 
temperatures in a system: T p , where formally the complete order parameter given by Eq.(|^) 
arises but its phase is a random quantity, i.e. < (p(x,y) >= 0, and another one, T BKT < T p , 
where the phase of the order parameter becomes ordered, so that < (p(x,y) >^ 0. It must 
be, however, stressed that the temperature T p is not a real critical temperature; it only 
denotes the characteristic region where the modulus of the order parameter achieves its 
maximal growth at T decreasing. Unlike T p the temperature Tbkt does correspond to phase 
transition when correlators < exp[i6(r, T)]exp[i6(r')] > as function of |r — r'| change their 
behaviour [ ffijfl . 

Let us find the n/-dependence of the temperatures T p and Tbkt- "Effective" temperature 
T p can be estimated from the set (p^6|)-(p!9|) in mean- field approximation by putting p = 0. 
Another temperature Tbkt follows from the Eqs. (17) and (19). 

It is impossible to solve the equations obtained analytically, so we shall do that by 
numerical calculation. Nevertheless, some assymptotical expressions for these temperatures 
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as functions of n/ can be found in analytical form: 

i) at ep/uJo — > one obtains Tbkt = £f/8, and T p satisfies the simple mean- field equation 
T p \ia(T p /eF) = u exp(— 2/A), where A = g^mjl^ is the dimensionless coupling constant. 

ii) ep — > 6jT; the critical point = e c f (or = nj, what is the same) at which 
T p = Tbkt = can be found from Eqs.(17)-(19). This unknown energy is the solution of 
the equation 



l _ A (W - e F )[(a - l)e F + u } _^ A (a - l)e F + w 
"2 n (a - l)ejp(W - e F + w ) ^ 2 11 (a - l)e F 



(20) 



(W is the conduction bandwidth determined by evident condition: W = k. 2 B /2m, where 
is the Brillouin wave vector). In other words, it follows from fl2!]| ) that because of the long 
waves phonons begin to quit the interaction transfer, there exist (at a > 1 only), the point 
where both T p and T BKT temperatures become zero what means that superconductivity (but 
not conductivity) is suppressed. Near this point this temperatures have the next behaviour 
(W -> oo): 



T — c cr 
J-BKT — 



1/5 



[a — 1 



144 v ' \ef[{a-l)e% + u ]J 
The results of numerical calculations of Eqs. (17)-(19) are shown on Fig.l. 
It is seen that due to the long- wave phonon suppression (in fact, to be supposed strong) 
the functions T p (nf) and Tbkt(^/) quickly on rif variation acquire a bell-like shape. This 
"bell" proves to be non-symmetrical, and its heigt, width and to some extent form depend 
on a(> 1). In some sense such a result is surprising because, as it is generally accepted, not 
long-, but short-wave intermediate bosons play the main role in attraction which appears 
in BCS-Eliashberg model due to electron-phonon (or any other boson) interaction. The 
sensitivity of superconducting properties of a 2D metal to long-wave part of the intermediate 
boson spectrum is rather unusual and allows to hope the more accurate consideration also 
results in a similar effect. 
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4 Conclusion 



The existence of two-different temperatures in underdoped HTSCs is now a well established 
fact. They, as it shown in many papers (see the rewiev [JTT|), are one of the consequences 
of two- dimensionality of their electronic and magnetic properties when the "ordering" of 
the order parameter modulus and phase takes place at different temperatures. In pure 2D 
system the temperature Tbkt, as was mentioned, has to be considered as critical one, and in 
the region T BKT < T < T p so called pseudogap, also normal, phase is formed in underdoped 
HTSCs. It is destroyed when T > T p or rif becomes so large (overdoped regime) that 
chemical potential of fermions and Fermi energy are indistinctive (/x = ep). And the aim of 
this paper was to demonstrate that in the superconducting system with indirect attraction 
the role of long-wave bosons (for instance, phonons, spin fluctuations and so on) can be 
crucial in the region where standard Fermi-liquid theory becomes already applicable. 

In spite of some qualitative similarity between experimental and obtained pictures the 
model considered is so simple (and even rough because the propagator (2) was, in fact, pos- 
tulated while it must contain the boson damping only) that any quantitative its use is almost 
senseless. Therefore it must be generalized on the taking into account such HTSC features 
as: quasi-two-dimensionality (which results in real T c ), intermediate boson dispersion uj(k) 
(as it takes place for spin fluctuations) and damping 7(k); pairing anisotropy, etc. These 
problems will be considered separately. 
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Figure 1: The characteristic patterns of the T — nf phase diagram of 2D metal with coupling 
constant A = 1. Doted and solid lines define the temperatures T p and Tbkt, respectively. 
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